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We present a theoretical study of lifetimes of interface states (IS) on metal-organic interfaces 
PTCDA/Ag(lll), NTCDA/Ag(lll), PFP/Ag(lll), and PTCDA/Ag(100), describing and explain¬ 
ing the recent experimental data. By means of unfolding the band structure of one of the interfaces 
under study onto the Ag(lll) Brillouin zone we demonstrate, that the Brillouin zone folding upon 
organic monolayer deposition plays a minor role in the phase space for electron decay, and hence 
weakly affects the resulting lifetimes. The presence of the unoccupied molecular states below the IS 
gives a small contribution to the IS decay rate mostly determined by the change of the phase space 
of bulk states upon the energy shift of the IS. The calculated lifetimes follow the experimentally 
observed trends. In particular, we explain the trend of the unusual increase of the IS lifetimes with 
rising temperature. 


I. INTRODUCTION 

Many well-defined interfaces between organic semicon¬ 
ductors and metals exhibit interface-specific electronic 
states that exist independently of the detailed molecule- 
substrate interaction^— Like the surface states (SSs) of 
clean metals they are a consequence of the breaking of 
translational symmetry perpendicular to the interface. 
In the Shockley-type interface states (ISs) that have been 
identified at Ag(lll) and Ag(100) interfaces the electrons 
are able to move almost freely parallel to the interface 
whereas the local charge density in the vicinity of the 
first molecular layer is strongly corrugated and resem¬ 
bles that of molecular orbitals. These general properties 
of the states have been revealed in recent experimental 
and theoretical works—For a few systems also the im¬ 
portant factors that determine their electronic structure 
could be investigated] 13 ! 16 

The dynamics of electronic decay and electron trans¬ 
fer processes at interfaces that involve these states, 
however, is not well understood. Previous two-photon 
photoemission (2PPE) experiments have measured life¬ 
times between 10 and 200 fs for electrons excited into 
normally unoccupied interface states above the Fermi 
level— ' 13 ’ 15 ! 16 From these short lifetimes a large overlap 
of the wave function with the metal has been concluded. 
These conclusions, although confirmed by density func¬ 
tional calculations ^ 10 ! 13 ~ — are based on very simplistic 
assumptions on the nature of electronic decay processes 
at such an interface. Many-body calculations, such that 
exist for surfaces states of clean metal surfaces^ have 
not been performed so far. 

In this publication, we make a first step in this di¬ 


rection. We perform a theoretical study of the decay 
of electrons in the interface states formed at the inter¬ 
faces of silver with ordered monolayers of such organic 
molecules as perylene-3,4,9,10-tetracarboxylic acid dian¬ 
hydride (PTCDA), naphthalene-1,4,5,8-tetracarboxylic 
acid dianhydride (NTCDA) and perfluoropentacene 
(PFP). First, on an equal footing we perform ab initio 
density functional calculations of the electronic structure 
of all the studied interfaces, of which the PFP/Ag(lll) 
interface is published for the first time here. Our calcula¬ 
tions show that the hybridization of molecular and metal¬ 
lic states is small in the region of the projected band gap 
of the metal. New elastic decay channels, which in princi¬ 
ple could open up due to reduced translational symmetry 
of the organic overlayers, will thus only have a weak in¬ 
fluence on the electron decay of the interface state. This 
allows us to focus on the inelastic decay of electrons ex¬ 
cited to the IS. We calculate the corresponding lifetimes 
in the self-energy formalism of many-body theory using 
the GW approximation!— In order to make the calcula¬ 
tions feasible, we use one-dimensional model potentials 
for an approximate description of the electronic struc¬ 
ture of the interfaces. The potentials are based on the 
so-called Chulkov potentia l 19 ’ 20 of clean surfaces. They 
are modified in order to reproduce the experimentally 
observed energy upshift of the Shockley surface state in 
the presence of the organic over layers. Comparison with 
the ab initio calculations allows us to judge how well 
the potentials reproduce the probability density of the 
states perpendicular to the interfaces. The results are 
in good overall agreement with the experimental data 
for PTCDA /Ag(lll), NTCDA/Ag(l 11), PFP/Ag(lll), 
and PTCDA/Ag(100). Calculated lifetimes are gener- 
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ally longer than experimental ones but agree well with 
experimental trends. 


II. AB INITIO CALCULATION OF 
ELECTRONIC STRUCTURE OF INTERFACES 


First, in order to have a detailed information about 
the electronic structure of the interfaces under study, 
we performed ab-initio calculations within the periodic 
slab geometry. We used the OPENMX (version 3.7) 
code^i, which is based on density functional theory and 
the linear combination of localized pseudoatomic or¬ 
bital (LCPAO) method.— We applied the general¬ 
ized gradient approximation (GGA) of Ref. [ 25 ] for the 
exchange-correlation functional. Also we exploited norm- 
conserving pseudopotentials^ in order to replace deep 
core potentials by shallow ones. For silver atoms, we 
set basis functions to Ag7.0— s2p2d2fl, while for hy¬ 
drogen, carbon, and oxygen atoms we use H6.0— s2pl, 
C6.0— s2p2dl and 06.0—s2p2dl, respectively. On the ex¬ 
ample of silver, this notation means that two primitive 
orbitals for each s, p , and d orbital and one primitive 
orbital for the / orbital were used for representation of 
the basis functions with the cutoff radius of 7.0 Bohr. 

To improve the description of the surface state (SS), 
we used the enlarged cutoff radius (9.0 Bohr) for silver 
atoms in the uppermost layers of the slab. The latter con¬ 
tains 10 silver layers together with the molecular mono- 
layer (ML) attached on one side of the silver film. The 
vertical distances between the Ag(lll) surface and the 
carbon atoms of NTCDA^ PTCDA^ and PFP— MLs 
are taken to be equal to its experimental values. The 
oxygen atoms of NTCDA and PTCDA monolayers on 
Ag(lll) are fixed at the same distance as carbon ones. 
The optimized geometry of PTCDA/Ag(100) found in 
Ref. [ij is used in our calculation. The real-space grid for 
numerical integration and solution of the Poisson equa¬ 
tion was specified by the energy cutoff of 250 Ry. The 
total-energy convergence was better than 0.027 meV. The 
surface Brillouin zone (SBZ) of the supercell was sampled 
with a 3 x 3 x 1 mesh of k-points. 

As seen in Fig. [I] due to the unit cell of the interfaces, 
which is larger than the (lxl) cell naturally used for 
the bare surfaces, the SBZ becomes smaller, and, conse¬ 
quently, the metal bands of the initial SBZ correspond¬ 
ing to the (lxl) unit cell get folded into the reduced 
SBZ, hybridizing at that with the orbitals of the molecu¬ 
lar ML. This leads to a surface band structure that does 
not exhibit the projected band gap at the F point any 
more. To restore the (1 x 1) unit-cell representation of 
the interface electronic structure, we performed an un¬ 
folding procedure using the BandUP code^lT— based on 
the method by Popescu and Zunger— Upon the unfold¬ 
ing procedure each electronic state mK (where m is the 
band index and K is the wavevector in the SBZ of the 
interface) is projected onto a set of corresponding points 
k i in the unfolded Brillouin zone (UBZ) of the (1 x 1) unit 
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FIG. 1: Band structure of (a) NTCDA@Ag(lll) in the re¬ 
laxed phase, (b) PTCDA@Ag(lll), (c) PTCDA@Ag(100), 
and (d) PFP@Ag(lll). The width of the fat bands reflects 
the extent of their localization in the vicinity of the molecu¬ 
lar ML. The red curves highlight the interface state originated 
from the Shockley surface state (resonance) of the respective 
bare silver surface. 


cell, resulting the weights W m K(k i ). The BandUP code 
deals with wavefunctions, expressed in the plane-wave 
basis set. Thus we perform calculations with the VASP 
code ; 34 ’ 35 based on the plane-wave basis, employing 
the projector-augmented wave (PAW) method. 36 ’ 37 The 
exchange-correlation was calculated within the GGA^ 
like in the LCPAO method. The energy cut-off was 
fixed at 350 eV and a 6x4x1 Monkhorst-Pack grid of 
k-points was used. These calculations are notably more 
time-consuming than the LCPAO-calculations, thus we 
restrict ourselves to the case of NTCDA@Ag(lll). 

The band-structure calculations indicate that the de¬ 
position of a ML on top of the silver surface influences 
the partly occupied SS in a qualitatively similar manner 
for different interfaces. In the cases of Ag(lll) surface 
the SS is shifted to a higher energy and transformed into 
the interface state (IS). In the case of the Ag(100) sur¬ 
face, the unoccupied Shockley resonance (SR) is also up- 
shifted in energy, which changes its character to a distinct 
electronic state. The magnitude of the upshift depends 
on the molecular type, coverage, surface orientation, and 
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adsorption geometr y 16 1 9 1 1 0 1 1 (see Table [TJ). Addition¬ 
ally, the bare-surface electronic structure is filled up by 
the molecular-derived weakly dispersive states, which are 
thought to have an effect on electron dynamics in higher 
lying states. 

In the unfolded bandstructure (Fig [2Ji) one can see 
a clearer picture how the deposition of an NTCDA ML 
modifies the electronic structure of the Ag(lll) surface. 
Below —3 eV one can see a set of d-bands, which are 
not affected by the molecular ML. At higher energies 
the bulk-derived s — p bands forming the gap at the T 
point at —1.2 eV < E < 3.6 eV, also have practically the 
same energies for surfaces with and without NTCDA. 
The main effect of NTCDA ML consists in the transfor¬ 
mation of the bonding and anti-bonding surface states 
(at -0.17 and +0.08 eV) of the bare 10-layer Ag(lll) slab 
into the SS of the clean side of the slab (at -0.06 eV) and 
the IS of the slab with NTCDA at Ejs = 0.5 eV. The 
latter value agrees with the result of LCPAO calculations 
Eis = 0.4 eV. The molecular states of NTCDA are un¬ 
folded to different points of the UBZ with small weights, 
thus forming a weak background. Note, that the spec¬ 
tral weight of the molecular-derived states is quite small, 
hence one could expect them to produce a rather minor 
contribution to decay processes. This question will be 
addressed in more detail later in Sec. m 


III. GW APPROXIMATION 

The lifetimes of ISs are calculated within the GW 
formalism^ As far as a fully ab initio calculation of 
lifetimes for the interface under consideration is a big 
challenge so far, one needs a model for the description of 
the electronic structure of these systems. Such a model 
can be based on the following propositions. First, the 
origin of the IS experimentally observed in the afore¬ 
mentioned interfaces is attributed to the upshifted SS 
(SR) of the bare surfaces. (See e.g. Refs SiEIIMl 
The properties of the resulting IS (the bulk penetration, 
near-surface localization, dispersion in k||, etc.) is sim¬ 
ilar to a Shockley-type state residing in the projected 
band gap of the bare (111) surface. Second, for IS the 
decay phase space is thought to be mainly provided by 
the projected Ag bulk states and is not affected by the 
band folding caused by using the large interface unit cell 
instead of the (lxl) unit cell of the bare Ag surfaces (see 
Fig. (2Jr). Third, due to the rather small spectral weight 
of the molecule-derived states in the unfolded electronic 
surface structure, their contribution to the decay can be 
neglected. We reduce thus the problem to a study of elec¬ 
tron lifetimes for a silver-like surface with a surface state 
modified in a way to reproduce the energy and dispersion 
of the considered IS. 

The mentioned simplifications allow us to use a one¬ 
dimensional (ID) potential mode l 19 ! 20 ’ 39 ^ 40 in description 
of the decay of IS electrons. In such a model the one- 
electron pseudopotential V ( z ) depends only on the z co¬ 


ordinate, being constant in the xy plane. The z axis is 
directed perpendicular to the surface outside the metal. 
The position of z = 0 corresponds to the plane of the 
topmost Ag atomic layer on the side where the ML is at¬ 
tached. The forms of the potential V(z) will be discussed 
further in sections UY10 

The method of calculation of lifetimes in the GW ap¬ 
proximation using a ID pseudopotential is described in 
detail elsewhere^ and here we give just a brief overview, 
indispensable for understanding of the further discussion 
of the results. Within this formalism the decay rate 
r e _ e of the electronic state in band i with wave function 
^ki(^, r ||) = +i(-) e * kr|1 , energy E ki = e+fc 2 /(2m*), and 
wavevector k is obtained as the projection of the imagi¬ 
nary part of the self-energy operator E onto this state: 


Fe-e = -2(4' k i|ImE|T ki ) (1) 

= — 2 ^^ J dzdz'M t j (z. z') 

3 

X j (27t) 2 t 1 ~ /<U'] — E qj) 

x ImW(z,z';k-q,E ki - E^). 


Here fqj is the Fermi factor and the wave-function prod¬ 
uct Mij(z,z') = ipi(z)tpi(z')ipj(z)tpj(z') is determined by 
the real eigenfunction ipi(z) being the solution of the one¬ 
dimensional Schrodinger equation 




ipi(z) = €ilfii(z) 


( 2 ) 


with the respective eigenvalue e;. In Eq. m, the self¬ 
energy is represented by the first term of the expansion 
in terms of the screened Coulomb interaction W, which 
is calculated within the random phase approximation. 
Thus, the many-body decay rate is determined by three 
main factors: (i) the phase space of the final states (qj), 
(ii) the overlap between the wave functions of the initial 
and final states, and (iii) the magnitude of the imaginary 
part of the screened Coulomb interaction ImW. The lat¬ 
ter is given in linear response theory by 


W(z, z'\ q, w) = v c (z,z'-,c\) + J dz 1 dz 2 v c (z,z 1 ;q) 

X x( z i, - 2 ; q, u)v c (z 2 , z'\ q), (3) 

where v c (z, z'\ q) = 2iTe~ q ^ z ~ z I /q is the 2D Fourier trans¬ 
form of the bare Coulomb interaction, and x is the 
density-density response function of interacting electrons, 
which is evaluated from the equation 


X(z,z'-q,uj) = x°(z,z';q,uj) + J dz 1 dz 2 X°( z , zi; q, w) 

x v c (zi,z 2 ;q)x(z2,z';q,v). (4) 

Here x°(ri ,r 2 ;oj) is the density-density response function 
of a non-interacting electron system: 

^k (fkj /k+qj ) Mij (z,z ) 


X 


] (z,z'; q,u) = 


(27r) 2 w + E kj - E k+ qi + irf 


( 5 ) 
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TABLE I: Experimental and theoretical values of the IS energy Eis (in eV) and lifetimes ns (in fs). Theoretical values of tis 
are given for the A-shifted / U-shifted scheme, accounting for inelastic electron-electron scattering only. 




PTCDA@Ag(lll) 

PTCDA@Ag(100) 

NTCDA@Ag(lll) 

PFP@Ag(lll) 

Eis 

Experiment 

0.57 ±0.02“ 

2.25 ±0.03" 

0.38 ±0.02“ 

0.1-0.2“ 


Theory 

0.55 

2.26 

0.40 

0.17 

ns 

Experiment 

53 ± 3“ 

3 < t < 18" 

115 ± 10“ 

— 


Theory 

110 / 270 

- / 24 

250 / 500 

1280 / 1850 


“Prom 

“Prom 

“From 


Ref. [H 


Ref. 
Ref. 3£j 


with r] being an infinitesimally small positive constant. 


IV. MODEL OF SHIFTED SURFACE STATE 

As noted above, the one-dimensional (ID) potential 
model may be used for the description of the decay of IS 
electrons. We start with the pseudopotential, introduced 
in Refs. 0 and 0 and a set of parameters that ensure 
the proper description of the surface electronic structure 
of the bare Ag surfaces. Further we need to modify the 
model to reproduce the energy of the IS, which is pre¬ 
sented by the SS, shifted to higher energies, together with 
the energy of the n = 1 IPS and the gap edges.. 

A way to achieve the shift of the SS, is to change the 
corresponding energies e*, entering Eqs. ([2]), ©. ([I]) and 
©” by hand”, while leaving the wavefunctions ipi ( z ) un¬ 
changed. Hereafter, we refer to such a scheme as to the 
one of the E-shifted surface state. In this case we do not 
change the overlap between the IS and the bulk states, 
while the phase space and to a certain extent the screened 
interaction are modified. Another way is to tune the pa¬ 
rameters of the ID potential in the near-surface region in 
a way, which provides a shift of the surface-state energy 
towards higher values at the f point (see Fig. [2jo,c). This 
scheme will be referred to as the one of the V-shifted sur¬ 
face state. Here, all the aforementioned factors that de¬ 
termine the inelastic decay are affected. In both schemes 
the energy of the IPS (n = 1) and the surface barrier are 
held unchanged. 

In Fig. [3l the starting-point ID potential and the cor¬ 
responding wave functions, as well as their counterparts 
at the largest energies considered, are shown. As seen in 
the figure, the shifted-SS penetration into the bulk be¬ 
comes smaller upon modifying the potential, leading to a 
redistribution of the SS charge density into the vacuum 
side. This redistribution to a certain extent reflects the 
situation observed in ab initio calculations, where in the 
case of the interface the charge outside the metal is larger 
than in the case of the bare surface (see Fig. a. 

On the one hand, the H-shifted scheme looks more con¬ 
sistent, because the energies et and wavefunctions 
are the solutions of the same Schrodinger equation ©. 
However, Fig. QJtop) shows the ab initio calculated wave 
function of the IS in comparison with the SS of the bare 
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FIG. 2: (a) Unfolded band structure of NTCDA@Ag(lll) 

(dots) in comparison with the bare Ag(lll) surface. The size 
of the symbol reflects the weight IF m K(k;) of the state in 
the unfolded BZ. The green (red) symbols denote the states 
with Wm K(k,) > 0.1 (W m K(ki) < 0.1). Sketch of the scheme 
realized in the model of the U-shifted surface state for Ag(lll) 
(b) and Ag(100) (c). The SS is shifted towards higher energies 
by modifying the ID pseudopotential of the bare Ag surfaces 
within the near-surface region. In the case of Ag(100), we 
start from the Shockley resonance (SR). 


silver surface, where along with the aforementioned redis¬ 
tribution to the vacuum side the penetration, and, con¬ 
sequently, overlapping with the bulk states, remains to 
be similar to the case of the bare surface. This means 
that the E-shifted SS scheme is valid and can be more 
suitable in the lifetime description. 

As an advantage, both schemes developed allow us to 
get F e _ e as a function of the IS energy, instead of specific 
calculations for a given IS energy. The decay rate of the 
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FIG. 3: One-dimensional potential and respective wave func¬ 
tions for the surface (interface) states and bulk states in the 
case of the Ag(lll) and Ag(100) surfaces at different energies 
of the shifted surface state ( E SS / IS ). The bulk charge density 
was normalised to fit the same scale as SS/IS. 


V-shifted SS electrons at f as a function of the SS(IS) 
energy is shown in Fig. [5](top) for the Ag(lll) surface. 
For example, for the experimental values of the IS energy 
observed at 300 K for NTCDA and PTCDA monolay¬ 
ers on Ag(lll) (see Table Q]) the decay rate F e _ e = 1.3 
meV (t ^0.5 ps) and 2.4 meV (r ~0.3 ps), respectively. 
Such a difference between the indicated values of T e _ e 
is caused by the fact that for the bigger IS energy there 
is a larger decay phase space formed by the bulk states. 
However, the difference could be bigger, if the aforemen¬ 
tioned overlap between the SS(IS) and the bulk states did 
not decrease with the increasing SS(IS) energy. Calcu¬ 
lations performed with the unmodified bare-surface ID 
potential but with changed SS(IS) energies (If-shifted 
surface state) clearly show how the decreasing overlap re¬ 
duces the gain in the decay phase space (see dashed line 
in Fig. EKtop)). In this case, for NTCDA and PTCDA 
monolayers on Ag(lll) the decay rates are T = 2.6 
meV (r ~0.3 ps) and 6.0 meV (r ~0.1 ps), respec¬ 
tively. The absolute values of the lifetimes found within 
both schemes for the interfaces are substantially longer 
than those experimentally observed: tI Jtg DA@Ag{m) = 



FIG. 4: Top: Ab-initio calculated wave functions of the in¬ 
terface and former LUMO states as compared with that of 
the surface state of the bare Ag(lll) surface in the case of 
NTCDA@Ag(lll). Bottom: The same wave functions as ob¬ 
tained within the ML-pseudopotential model (see the text). 


115 ± 10 fs and Tpf’<y£> J 4 @A ff (iii) = 53 ± 3 fs at 300 K and 

T N S TCDA@A g (iu) = ^3 fs and Tpp CDA@Ag ^ 11Vj = 26 fs at 
90 K as reported in Ref. |Tj. However, their ratio is very 
close to its experimental counterpart. Moreover, the pre¬ 
sented results allow us to nicely reproduce the trend of 
the unusual increase of the lifetimes with rising temper¬ 
ature (see Fig. Otop)). Actually, at higher temperature 
the absorption distance becomes larger and, as a con¬ 
sequence, due to the weakened interaction between the 
molecular ML and metal substrate, the IS energy gets 
smaller, providing by this the smaller decay rate. 

For the PFP@Ag(lll) interface no experimental data 
is available on the IS lifetime. Our calculations yield a 
very small inelastic decay rate ~ 0.5 meV, corresponding 
to the lifetime of 1-2 ps due to the low energy of the 
IS, and hence small phase space for the inelastic decay. 
However, this value should be sufficiently shortened by 
elastic decay channels, e.g. electron-defect scattering. 

In the case of the Ag(100) surface, the decay rate of 
the surface resonance (SR) cannot be calculated within 
the present framework. However, when the SR is pushed 
into the band gap, it becomes the IS. The decay rate of 
the latter decreases with its energy (see Fig. [^bottom)) 
in spite of the growing decay phase space. Such a de¬ 
crease is caused by the rapid reduction of the IS-bulk 
states overlap (see Fig. n. For the IS energy that corre¬ 
sponds to the PTCDA/Ag(100) interface, the decay rate 
is of 27 meV. The corresponding lifetime r = 24 fs is 
surprisingly close to the upper limit of 18 fs found ex¬ 
perimentally in Ref. [l5]. Note that we consider here the 
R-shifted SS scheme only, since we cannot extrapolate 
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FIG. 5: Calculated decay rate T of the shifted surface state 
as a function of its energy. Solid (dashed) lines correspond 
to the V-(.E-)shifted scheme, respectively. Light-gray stripes 
cover the energy intervals, which correspond to the experi¬ 
mental values of the IS energy, including its variations with 
temperature. In the case of Ag(lll), red points show the de¬ 
cay as obtained within the E-shifted scheme with taking into 
account the contribution of the former LUMO. Error bars 
here are caused by changing the former LUMO energy (see 
Fig- 0 within the ±50 meV interval. 


the resonance-like wavefunction of the bare Ag(100) sur¬ 
face to the projected band gap energy region. 


V. MONOLAYER PSEUDOPOTENTIAL 

In an attempt to take into account a possible contribu¬ 
tion coining from the molecular-derived states presented 
in Fig. |T]by the fat bands (or red circles in Fig. the 
localization of these states lying energetically below the 
IS (the former LUMO in the Ag(lll) interfaces) should 
be reproduced. (See Fig. |D(top)) For this purpose we de¬ 
velop a pseudopotential that models the molecular ML 


as a quantum well and contains a barrier, separating the 
ML from the bulk Ag(lll). We refer to this model as 
ML potential. This model allows us to get the shifted SS 
and the former LUMO with quite close energies (repro¬ 
ducing the ab initio values) and localized in the interface 
region and at the molecular layer, respectively, with the 
overlapping f \<pis{z)\ 2 \ipLUMo{z)\ 2 dz as obtained from 
the ab initio calculations. 

In Fig. HKbottom), we plot the obtained model wave 
functions. Note, that in the calculation of lifetimes we 
have replaced the wavefunction of LUMO </?lumo (z) by 
Plumo (z) / y/N in order to take into account the unfold¬ 
ing onto the (1 x 1) Ag(lll) BZ. Here N is the num¬ 
ber of atoms in the Ag(lll) surface layer unit cell of 
the interface: N = 33 and 24 for PTCDA and NTCDA, 
respectively. The resulting contribution accounting for 
transitions from the SS(IS) to the LUMO was estimated 
to be of ~ 0.3 — 0.7 meV. Fig. [Sjtop) shows the quite 
moderate effect of this contribution being added to the 
decay rate of the ^-shifted SS scheme by red dots. By 
error bars, we demonstrate how variations in the former 
LUMO energy (±50 meV) can affect the resulting T. 

As one can see, the calculated lifetimes are generally 
longer, than the experimental ones. On the one hand, 
this is a typical picture, because there are usually some 
other decay channels, which are not taken into account 
in the model. On the other hand, the difference between 
theoretical and experimental results may be caused by 
inaccurate description of the charge density distribution. 
As follows from ab-initio calculations, in the interfaces 
under study the presence of the molecular ML causes a 
redistribution of not only the SS charge density in the 
vicinity of the ML plane, as we show above, but also 
the bulk-states charge density. This simultaneous redis¬ 
tribution can produce an additional overlap between the 
IS and the bulk states, making the IS lifetime shorter. 
In order to take it into account, one should modify the 
ID pseudo-potential in a way to properly describe its 
behavior within the ML region. Since it involves the 
bulk states and affects the image-potential tail of the ID 
pseudo-potential, it should be done consistently with a 
study of energies and lifetimes of image-potential states, 
which, in turn, need to be experimentally analyzed first. 
On the other side further experimental studies of IS life¬ 
times for the organic monolayers on top of the Ag(lll) 
surface might shed light on peculiarities of this surface, 
that provide a distinct decay process as compared to the 
Ag(110) surface. 


VI. CONCLUSIONS 

We have presented a theoretical study of lifetimes of 
interface states (IS) on metal-organic interfaces within 
one-dimensional pseudopotentials. The presented results 
allowed us to address the question about the description 
and explanation of recent experimental data. We have 
demonstrated, that the folding of the BZ of the Ag sur- 
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face due to sustantial enlarging of the surface unit cell 
upon the deposition of organic molecules does not have 
a drastic effect on the phase space of final states for elec¬ 
tron decay. 

By means of the unfolding procedure we demonstrated, 
that the BZ folding upon organic monolayer deposition 
plays a minor role in the phase space for electron decay. 
Actually, it introduces only a weak background to the 
bandstructure of the (1 x 1) unit cell, and hence weakly 
affects the lifetimes. The presence of the unoccupied 
molecular state below the IS gives a small contribution to 
the decay rate of the IS, while the lifetime is mostly de¬ 
termined by the change of the phase space of final states 
upon the energy shift of the IS. 

In the case of PTCDA@Ag(100), the IS lifetime ob¬ 
tained in the E-shifted SS scheme agrees well with the 
experimental data, while for PTCDA and NTCDA at 
Ag(lll) these model strongly overestimates the corre¬ 
sponding lifetimes of the IS. Being applied to these 
Ag(lll) interfaces, the E-shifted SS scheme, which is 
based on the phase-space description, yields shorter life¬ 


times, but still quite long as compared with the ex¬ 
periment. However, our calculations provide the ratio 
T N S TCDA@Ag(ui)/ T PTCDA@Aa(iii) in agreement with the 
experiments, and explain the trend of the unusual in¬ 
crease of of the IS lifetimes with rising temperature. 
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